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Abstract 



We find an analytical expression for the specific heat of a Fermi gas in a 
harmonic trap using a semi-classical approximation. Our approximation is 
valid for fcgT > hw X) y^ z and in this range it is shown to be highly accurate. 
We comment on the semi-classical approximation, presenting an explanation 
for this high accuracy. 
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Trapped ultra-cold Fermi gases have received a lot of attention in the last few years. 
This was partly triggered by the first observations of Bose-Einstein condensation in 1995 
|jl|-0|. Since then, the field has vastly expanded. The trapping and cooling of fermions is 
in a much less developed state than that of bosons but it will probably not be long before 
quantum degeneracy is achieved [|J. (For recent experiments see |5]||.) 

On the theoretical side, Butts and Rokhsar have calculated the specific heat nu- 
merically and provided analytical results for the spatial and momentum distribution of a 
non-self-interacting spin polarised Fermi gas at T = in a harmonic trap using the semi- 
classical (Thomas-Fermi) approximation. Their results are valid for large particle numbers. 
In addition, Schneider and Wallis || made a numerical study of a similar gas and focused 
on the effects of small particle numbers. The present paper complements that of Butts 
and Rokhsar in that we provide an analytical expression for the specific heat of a non-self- 
interacting spin polarised Fermi gas in the context of a semi-classical approximation. We 
then compare it to the exact result given by a numerical calculation. As we shall see, the 
semi-classical approximation produces extremely good results for kT > hw Xt y )Z , where w x>y ^ 
are the frequencies of the trap in the three spatial directions, for particle numbers as low 
as N = 1000. This high level of accuracy is perhaps unexpected. We will discuss this and 
an explanation will be given. The characteristic temperature for quantum degeneracy, the 
Fermi temperature, is given by kTp = hw(6N)3 for an isotropic trap. The condition we 
impose is therefore valid for a quantum degenerate regime. It is also worth noting that this 
condition holds in current experiments with trapped fermions ||. 

The effects of interactions has been dealt with in several works [PHT^ and in particular 
the possibility of a BCS transition has been studied in the case of 6 Li 0-11 . Considering 



the gas to be non-interacting is a very good approximation for a dilute neutral atomic gas. 
The only case where interactions could be important is trapped 6 Li with two hyper-fine spin 
states fl2|]. 6 Li has an anomalously large s-wave scattering length. At least two different 



states are needed for s-wave scattering as it is forbidden for particles in the same state. 



However, the results of ref. [TT[] are for kT << hw, lying outside the range of temperatures 
concerned in the present work. 

Consider a gas of fermions in a harmonic potential. Neglecting interactions between the 
particles, the energy levels of each particle are 

E nin2 n 3 = (ni + ^)ftu>i + (n>2 + 7p^2 + (n 3 + ^)fiw 3 , (1) 
rii = 0, 1,2, . . . , 

where u>i,2,3 are the frequencies of the trap. 

We will use grand-canonical statistics throughout. This makes the calculations much 
easier and is justified as the difference between canonical and grand-canonical results is 
minute for the particle numbers we consider ||. 

The number of particles in the system is 

N = Y\e^ En -ti + I]" 1 . (2) 

n 

The internal energy of the system is given by 
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U = Y,E n [e^ E ^ + . (3) 

n 

H is the chemical potential and the sum in (2) and (3) is over all particle states. Here it is 
convenient for simplicity of the expression, to make the substitutions 



Wi 


= w , 


w 2 


= Xw , 


w 3 


= X'w , 


X 


= f3hw , 


fi 


/ 1 

— hw I- 



~ 2 «J , (4) 

where x, e, A and A' are newly defined dimensionless variables. The internal energy becomes 

1 1 
U = -hw(l + A + X')N + -hwu , (5) 

where 



oo oo oo 



« = 2EEE k[e^ + I]" 1 , (6) 

ni=0 ri2=0 n3=0 

k — rt\ + An 2 + A'n3 . 

Only the second term contributes to the specific heat since iV is held fixed. We have 

_ (dU\ 1 (du\ 

\ / N,w \ / Ai,u> 

Note that the partial derivative at constant particle number and trap frequency is the natural 
way of defining the specific heat for a gas in the situation concerned. After a change of 
variables from T to x this becomes 

where ks is Boltzmann's constant. 

The simplest approximation consists in replacing the triple sums in (6) by a triple inte- 
gral. Thus, 

r+oo r+oo r+oc 

u~v = 2 / / k[e x ^ + l]- 1 dy 1 dy 2 dy 3 (9) 

JO JO Jo 

and k is now the real function k — y 1 + Xy 2 + X'y 3 . Changing the integrating variable to k 
we have 



1 f+oo 

V= XX ' 



r+oo 

k 3 [ e x ^ + l]- 1 dk . (10) 

Jo 



Note that for the isotropic case we could have transformed the triple sum in (6) into a 
single sum and only then convert it to an integral, yielding 
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1 f + OO 

v = — (k 3 + 3k 2 + 2k)\e x{k+e) + l}~ 1 dk , (11) 
AA' Jo 

which, at first sight, might be thought to be a more accurate expression. For the anisotropic 
case, though not trivial, we could likewise have a density of states expanded to three terms 
17] , 18 ]. We will coment on this shortly. 

To solve the integral in (10) we Taylor expand the integrand in powers of e x<yk+e ^ and 
integrate each resulting term. Here we have to consider two different cases: e > and e < 0. 
For e > this yields 

a oo —nxe 

7-iE(-D" +1 V- < 12 ) 



From (8) and (12) and retaining the fact that v is an approximated version of u, we have 

C 12 , 3 fd(ex)\ 

S"^^ + AAv(Vj N ,/ 3 ' (13) 
where A^ denotes the sum 

oo — n:r|e| 

^ = E(-i) n+1 — -• (14) 

n=l U 

The use of the absolute value of e is redundant in the case we are considering (e > 0) but it 
is useful when we consider the e < case. 

To obtain an expression for the partial derivative in (13) we use the equality 



ox J 



= , (15) 

N,w 



with N, the number of particles, being 



N~-4h- (16) 

AA'x 3 v ; 

Expression (16) is obtained using the same process as was used to obtain (13). We then 
have 



and 



( 3 ^4-3^1 • (18) 



k B w x 3 v M, 

Finally we note that the sums Ai can be put in terms of polylogarithms [BJ, yielding 



AA'— ~ 3x" 3 
k B 



l r , . JLi 3 (e-™) - lLi 3 (e- 2 ™)) 2 



4Li 4 (e- xe ) - -Lu(e- 2xt ) - 3 



2 4V ' Li 2 (e-**) - \Li 2 {e~ 2 ^) 



(19) 
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For e < 0, the procedure is similar but slightly more complicated. The integral of 
expression (10) must be divided in two parts, according to whether the exponential is greater 
or smaller than 1, as the Taylor expansion of the integrand is different for the two cases. 
The final result for e < is 

AA'f ~ -12A 4 x" 3 + 21C(4)x- 3 + 6C(2)e 2 x- 1 + Ux + {A \~ ^ "sff f , (20) 

where ( is the Riemann zeta function. Putting the sums Ai in terms of polylogarithms we 
have 

AA'— ~ [-12Lt A (e xe ) + -Li A (e 2xe ) + 21C(4)]aT 3 + 6C(2)e 2 x~ 1 + -e 4 x + 
kb 2 2 

[Li 3 (e") - jLi 3 (e**) - ((2)ex - |e 3 a; 3 ] 2 

L^(e-)-|Li 2 (e 2 -)-C(2)-|e 2 a; 2 ' 1 ' 

Note that these analytical expressions are not an approximation over the Thomas-Fermi 
result. They are the exact Thomas- Fermi result, which of course, is in itself an approxi- 
mation. We computed the specific heat numerically using an exact expression composed of 
several sums, which is easily obtained inserting (6) in (8) and using (15) to get an expres- 
sion for (d(ex)/dx)N,w This is then compared to the values for the specific heat given by 
(18) and (20). The accuracy of the Thomas- Fermi approximation is extremely good for x 
of order 1 or less, as can be seen in figures 1 and 2. As an example, for N = 1000 and 
X\ = X2 = £3 < 2.5 the error in C is always less than 1% and if we take Xj < 1, the error is 
always less than 0.1%. For larger values of x the approximation rapidly deteriorates. The 
results for a higher number of particles are even more accurate, as expected from this kind 
of approximation [[/]. 

Given that for this approximation we have replaced triple sums by triple integrals, such 
a level of accuracy is quite surprising. Note that the expression inside the sums in (6) falls 
off exponentially. Therefore the first few terms are very important and expression (9) should 
differ significantly from (6). And indeed it does. The reason our results are so accurate is 
that the value of e used in (18) and (20) is not the same as the one used in (6). For our 
Thomas- Fermi approximation we took e from expression (16), ie, an approximate e. So, the 
insertion of an approximate e in an approximate expression leads to a fortunate cancelation 
of errors. We have tested this by inserting the exact value of e into (18) and (20). This 
leads to much worse results. Also, we have tried to further improve our approximation by 
considering the first two terms of the cubic polynomial in (11) instead of only the first term 
and doing the same to the number of particles, obtaining an additional term in (16) as well 
as additional terms in (18) and (20). To our surprise, this only increases the error. We can 
only assume that the very high level of accuracy of the Thomas-Fermi approximation is a 
happy coincidence. In fact, if we take the exact value of e and input it in expressions (18) 
and (20) and in the equivalent expressions obtained using k 3 + 3k 2 in (11) instead of only 
k 3 , then the latter yield the most accurate result, though still not a very good one. 
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FIGURES 



FIG. 1. the familiar curve of the specific heat ploted against x = Hw/ksT for A = A' = 1 
(isotropic harmonic oscilator) and N = 1000. Both the analytical and numerical results are shown, 
but they are so similar that it looks like one curve only. 

FIG. 2. the ratio C 'analytical i 'Cnumericai ploted against x for A = A' = 1 and N = 1000 
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